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ABSTRACT 

X-ray emission from many rotation-powered millisecond pulsars (MSPs) is observed to be of predom- 
inantly thermal nature. In PSR J0437-4715, the nearest MSP known, an additional faint power-law 
tail is observed above 2.5 keV, commonly attributed to non-thermal magnetospheric radiation. We pro- 
pose that the hard emission in this and other similar MSPs is instead due to weak Comptonization of 
the thermal (blackbody or hydrogen atmosphere) polar cap emission by energetic electrons/positrons 
of small optical depth in the pulsar magnetosphere. This spectral model implies that all soft X-rays 
are of purely thermal origin, which has profound implications in the study of neutron star structure 
and fundamental pulsar physics. 

Subject headings: pulsars: general — pulsars: individual (PSR J0437-4715) — 
stars: neutron — X-rays: stars 
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As the nearest and brightest rotatio n-powered mil- 
liseco n d pulsar (MSP), P S R J04 37-4715 l|,Tohnston et all 
119931: Ivan Straten et all 1200 Ij) has been the best 
studied in the X-ray domain. ROSAT, Chandra, 
and XMM-Newton observations 



l|Zavlin & Pa vlov 199S, 

51 120021: lZavhnll200fi^ have revealed that the 
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emission in the 0.1-10 keV band consists of two thermal 
components, and a faint power-law (PL) tail, with pho- 
ton index F 2. The latter can only be discerned in the 
spectrum beyond 2.5 keV. Due to the limited spectral 
information above ~3 keV, the nature of this hard X-ray 
component is unclear. 

In this Letter, we examine the plausibility of various 
emission mechanisms capable of producing the hard X- 
rays observed from J0437-4715 (henceforth J0437). The 
paper is organized as follows. In §2 we describe the ob- 
servations used and their analysis while §3 outlines vari- 
ous interpretations of the emission from J0437. In §4 we 
present the results of our analysis and offer a discussion 
and conclusions in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

We have retrieved archival Chandra and XMM-Newton 
observations of J0437 to investigate the X-ray emission 
from this MSP. The Chandra ACIS-S observations were 
performed on 20 00 May 29-30 (25.7 ks effective exposure, 
3.9' off axis; see iZavlin et al.ll200l for details). We first 
reprocessed the level 1 event files using the most recent 
calibration release and the CIAO 3.3 software package^. 
The data were then filtered to the nominal Chandra en- 
ergy band (0.3-8 keV). The source counts were obtained 
from a 5" circular region centered on the radio position of 
the MSP. This circle encloses >90% of the total energy. 

The XMM-Newton data consist of a 70-ks observation 
with the EPIC-MOSl /2 detectors in full-frame mode and 
EPIC-pn instrument in fast timing mode obtained on 
2002 October 9 (see [ZavUn 2006). The raw data were 
reprocessed using the emchain and epchain pipelines of 

^ Chandra Interactive Analysis of Observations, available at 
[htt p: / /asc. harvard.edu/ciao/ 



the SAS 6.5.0 software^ and were screened for instances 
of high background flares, which were subsequently re- 
moved. In the case of MOSl/2, the source counts were 
extracted from 70" circles, which enclose >90% of the to- 
tal energy. To minimize background contamination the 
MOSl /2 data were restricted to the 0.3-8 keV band. For 
the pn detector, the MSP counts were extracted from 
CCD columns 28-39 and were fihered to 0.5-10 keV, 
where the detector response is best known. 

For all observations, the background was obtained from 
a source-free region near the position of the MSP. The 
spectral fits were performed in XSPEC 11.3.1'^. 

3. ORIGIN OF THE HARD X-RAY TAIL 

The inherent faintness of the hard PL component from 
J0437 and the limited bandwidth over which it can be 
observed pose an obstacle in deciphering the origin of 
this emission. Nonetheless, using the available data we 
can by gain more insight by examining the plausibility of 
the various mechanisms for production of the observed 
hard X-ray tail. 

3.1. Non-thermal Emission 

The most widely accepted model for production of a 
PL spectrum in pulsars involves non-thermal radiation 
by relativistic electrons and positron s (e^) iri the p ulsar 
magne t osphe re . The analyses by Za vlin et alJ 1)2002') and 
iZavlinI l)2006j) of Chandra and XMM-Newton observa- 
tions found that the hard X-ray tail seen in J0437 is well 
described by a PL with spectral index F = 2.0 ±0.4. Un- 
der the assumption that this emission is of non-thermal 
magnetospheric origin, the PL is expected to extend from 
the optical/UV and beyond the hard X-ray range (>10 
keV). However, Hubble Space Telescope FUV observa- 
tions have revealed that extrapolation of the powerlaw 
to lower energies is inconsistend with the FUV flux from 
the MSP, with the lower bound F » 1.6 only marginall y 
consistent with the FUV flux l|Kargaltsev et al.ll2004j) . 
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The measured FUV flux is otherwise consistent with that 
of ^10^ K thermal emission from the whole NS surface. 
The apparent discrepancy casts doubt on the validity of 
t he PL model. 

iBecker fc Triimpeil (11999*) have proposed that the 
X-ray emission from J0437 and other nearby MSPs 
(e.g. PSRs J0030+0451 and J2124-3358) could be due 
to non-thermal processes in the magnetosphere. How- 
ever, based on more recent observations of these MSPs, 
we conclude that this interpretation is untenable. If the 
emission were indeed of purely non-thermal character, 
then such a broken PL spectru m grossly overestimates 
the obser ved optic al /UV fl uxes llKoptsevich et aLlEoO^t 
iMignani fe'Eecked 1:2003: K argaltsev et alJ l2004|) unless 
yet another break exist below 0.1 keV. No such behav- 
ior is observed in the X-ray spectra of other pulsars, in 
particular the luminous MSPs B1937-H12 and B1821-24, 
which show purely non-thermal pu lsed radiation (see e.g., 
IBecker et al]l200a IMineo et al.l l2004') . 

An alternative source of non-thermal X-rays is nebular 
synchrotron emission from an intrabinary shock formed 
due to interaction between the wind from the MSP and 
the companion star. Such a shock is believed to be 
present in other binary MSPs, namely PSRs B1957-I-20 
(iStappers et alJ l200^. J0024-7204W ijBogdanov et all 
|20nir"and~jT740-5340 (Grindlay et al. in prepara- 
tion). However, for J0437, the strength of the pulsar 
wind E 4: X 10^'^ ergs s~^ and the orbital separation'' 
of a = 1.2 X 10^^ cm do not favor this interpretation as 
the expected X-ray luminosity from the resulting shock 
is substantially lower than observed. Specifically, assum- 
ing a radius ~1.5 x 10^ cm for the 0.24 M p, white dwarf 
(WD) companion ijvan Straten et al]l200lD . we find that 
the total incident power on the face of the WD from 
the pulsar wind is ■^lO^^ ergs s~^, more than two or- 
ders of magnitude lower than the observed luminosity 
of the hard emission {Lx « 7 x 10^® ergs s~^ for 0.1-7 
keV). The large discrepancy can be reduced in a more 
contrived scenario, involving a highly anisotropic wind, 
confined to the orbital plane in a beam covering <0.15% 
of the sky. Alternatively, the companion may be bloated 
or surrounded by a large tenuous envelope of gas that 
interacts with the MSP wind. However, this assertion is 
in conflict with the observed optical properties of J0437, 
which are entirely consistent with that of a ~0.2 M© He- 
WD (see e.g.,[Bailvn 1993). Therefore, we conclude that 
an intrabinary shock is not a significant source of X-ray 
emission in the J0437 system. 

3.2. Thermal Emission 

In a purely thermal scenario the hard emission from 
J0437 may originate from a very small hot spot on the 
polar cap with a radius of order a few meters, radiating at 
temperatures in excess of 8 x 10® K. However, no current 
MSP polar cap heating model predicts such peculiarly 
high temperatures and small emission areas, making this 
model rather implausible. The validity of this model can 
be tested by deep observations at >3 keV. 

3.3. Comptonized Thermal Emission 

" In lZavlin etan l200l) the semi-major axis of the MSP (ai) is 
erroneously quoted as being the distance between the MSP and its 
companion (a = ai + 02). 
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Fig. 1.— {top) Chandra ACIS-S and XMM-Newton EPIC- 
MOSl/2 and EPIC-pn spectra of J0437 fitted with a two temper- 
ature Comptonized blackbody model (compbb). The lower panel 
shows the best fit residuals, {bottom) Model spectra used in the 
fit. The solid line is the total spectrum, while the dashed and dot- 
dashed lines show the individual emission components (see text for 
best fit parameters). 



The hard X-ray tail can be naturally produced by re- 
processing of the thermal emission via inverse Compton 
scattering (ICS) by an optically thin (r < 1) hot plasma. 
As a thermal photon emitted at the NS surface prop- 
agates through the tenuous co-rotating plasma in the 
pulsar magnetosphere, it may undergo ICS by energetic 
e^. In this regime, repeated scatterings of the thermal 
seed photons result i n a P L distribution of photo ns (see 
iRvbicki fc LightmanllT979t iNishimura et alJll98^ . Due 
to the low r of the magnetospheric plasma, only a small 
fraction of the thermal photons are scattered, resulting 
in a faint Comptonized tail. 

A major advantage of this model is that, unlike the 
case of a non-thermal (magnetospheric) PL, it does not 
extend below ~1 keV (see Fig. 1), thus, aleviating any 
discrepancies with optical/UV data. Therefore, based 
on the available observational evidence and the relative 
simplicity of the Comptonization model, we deem this 
interpretation the most plausible. 

4. RESULTS 

We have examined the effect of inverse Compton scat- 
tering (ICS) on the thermal emission from MSPs for both 
a blackbody (BB) and neutron star hydrogen atmosphere 
(NSA, see Roniani 1987; Zavhn ct al. 1996) model. For 
the BB case we employed the compbb model in XSPEC 
(Nishimura et al. 1986). This particular model assumes a 
plane-parallel, semi-infinite scattering medium with the 
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Fig. 2. — (top) Example Comptonized BB model spectrum for 
fcTeff = 0.23 keV, fcTe = 100 keV, and r = 0.1. The dashed line 
shows the initial thermal photon distribution. The inset shows 
example 1-, 2-, and 3cr contours for r versus kTe to illustrate the 
correlation between the two parameters, (bottom) Comptonized 
NSA model IHeinke et al.ll2(K!3) spectrum for T^s = 1.3 X 10*^ K, 
kTe = 100 keV, and t = 0.1. The photon energies have been 
corrected for the gravitational redshift assuming i? = 10 km and 
M = 1.4 M0. The dotted lines correspond to different viewing 
angles relative to the surface normal (cosS = 0.5, 0.25, and 0.1, 
from top to bottom, respectively), while the solid line is for cos 6 = 
1. In both plots the dashed line shows the initial thermal photon 
distribution. The flux normalization is arbitrary. 



source of thermal photons at the bottom. It is valid 
over the range E <§; kT^, < mc^, where E is the ini- 
tial photon energy, is the scattering e^ temperature, 
and mc^ is the electron rest energy. Figure 1 shows the 
X-ray spectrum of J0437 fitted with a two-temperature 
Comptonized BB model using compbb. Throughout the 
spectral fits we fixed the hydrogen column density to- 
wards the MSP to iVn = 2 X 10^^ cm"^ and electron 
temperature to — 150 keV. The best fit parameters 
we obtain are Ti = (2.93±0.03) x 10^ K, i?i =44±llm, 
Ti = 0.06±0.02, T2 = (I.16±0.01) X 10^ K, R2 = 300±40 
m, and T2 = 0.09±0.0I. The corresponding luminosity in 
the 0.3-10 keV is found to be Lx « 3 x 10^" ergs s^^. AU 
uncertainties quoted above are Icr. The minimum value 
of xt = 1-36 (for 464 degrees of freedom) was found to be 
identical to that in the case of the BB-I-BB-I-PL model. 
Given the uncertainties in the cross-calibration of the dif- 
ferent detectors, this value of xt indicates an acceptable 
fit. As is generally the case for ICS, Tg and r are strongly 
correlated (see inset in Fig. 2) so a good fit was obtained 
for a wide range of kT^ and r. 



The compbb model assumes a thermal distribution of 
scattering e''^. In reality, these particles probably follow 
a non-thermal (PL) distribution over a certain range of 
energies. Notwithstanding, the exact energy distribution 
of the e^ has no bearing on the validity of this model 
since in both cases, multiple scatterings result in a PL 
distribution of photons. 

To study the effect of Comptonization on the emis- 
sion from a NSA we used the unmagnetized hydrogen 
atmosphere model from Heinke et al. (2006) (see also 
McClintock et al. 2004) and the Comptonization algo- 
rithm of iNishimura et al. (1986). As expected, we have 
found that the NSA model yields a qualitatively similar 
spectrum to that of the BB case (see Fig. 2). However, 
crucial differences do exist, steming from the inherently 
anisot ropic and energy d ependent emission pattern of a 
NSA llZavlin et al.lll99l . In particular, at a given view- 
ing angle, the intensity of the NSA radiation is subject to 
limb darkening as well as a projection effect. The scat- 
tered photons, on the other hand, probably have a differ- 
ent angular distribution. This suggests that depending 
on the viewing geometry, the observed flux may show 
variations in the relative contribution from the Comp- 
tonized component. 

It is important to emphasize that the exact geome- 
try and location of the scattering medium relative to 
the MSP polar caps is unknown and cannot be deter- 
mined from the time integrated spectra alone. In prin- 
ciple, these properties can be inferred by studying the 
modulation of the Comptonized radiation as a function 
of the rotational phase of the MSP. To investigate the 
temporal behavior of the PL tail, we used the timing 
data from the XMM-Newton EPIC-pn detector. Figure 
3 shows lightcurves of J0437 folded at the MSP spin pe- 
riod {P — 5.76 ms) for the 0.3-1 keV and 3-6 keV energy 
ranges. In the Comptonization scenario, these two bands 
contain purely (unscattered) thermal and just Comp- 
tonized photons, respectively. Unfortunately, due to the 
limited count statistics and the high instrument and sky 
backgrounds, which dominate the emission for >3 keV, 
at present, we cannot provide any meaningful constraints 
on the properties of the Comptonizing medium. We are, 
however, able to place a limit on the pulsed fraction of 
<50% for >3 keV. This value further constrains a non- 
thermal magnetospheric origin since such emission is ex- 
pected to be almost entirely pulsed. 

5. DISCUSSION AND CONCLUSION 

We have demonstrated that the X-ray emission 
from the nearby MSP J0437 in the 0.1-10 keV band 
is well described by a two-temperature Comptonized 
thermal spectrum. This model provides an attrac- 
tive and relatively simple alternative to the non- 
thermal magnetospheric em ission model. Since in the 
'Harding fc Muslimovl l|2002D model of high -energy emis- 
sion from MSPs, ICS is the primary mechanism for pair 
production and polar cap heating, the presence of a hard 
Comptonized tail would provide a strong confirmation 
that this mechanism does indeed operate in MSP mag- 
netospheres. 

We note that the shape of the spectra in the non- 
thermal magnetospheric and Comptonization cases are 
indistinguishable in the 0.1-10 keV band. The primary 
difference arises in the optical/UV range, where only 
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Fig. 3.— XMM-Newton EPIC-pn lightcurves of J0437 in the 0.3- 
1 keV (top) and 3-6 keV (bottom) bands. The dotted line in the 
bottom panel shows the background level. Two cycles are shown 
for clarity. 



as PSRs J0030-I-0451 and J2124-3358, would provide a 
more definitive test of the two interpretations (the exis- 
tence of a hard X-ray component cannot be ascertained 
in present X-ray observations of these MSPs due to the 
Hmited counts). For these soHtary MSPs, the absence of 
a stellar companion has allowed ultra-deep optical obser- 
vations. which have found no emission down to V ^ 28 
ijKoptsevich et al.ll2003t iMignani fc Beckeill2003|) . 

The Comptonization model proposed in this Letter im- 
plies that all of the soft emission (<1 keV) from J0437- 
4715 is purely thermal (unscattered) radiation. This has 
profound ramifications in the study of fundamental NS 
properties since thermal emission from the surface of neu- 
tron stars reveals important information regarding the 
basic properties of this class of compact objects. Pavlov 
& Zavlin (1997) have shown that X-ray spectral and tim- 
ing observations of MSPs can be used to measure funda- 
mental NS parameters such as the mass-to-radius ratio 
[M/R) of the underlying NS. In addition, the polar caps 
represent "footprints" of the pulsar magnetic field, al- 
lowing constraints on the global field geometry. Finally, 
since the properties of the hard emission component are 
determined by the energy distribution of scattering par- 
ticles and their optical depth, study of this radiation can 
also serve as a probe of the particles populating the pul- 
sar magnetosphere. An attempt to constrain the proper- 
ties of J0437 outlined above is currently under way and 
will be presented in a future paper. 



the former is expected. As mentioned previously, for 
J0437 the non-thermal model has been found to be only 
marginally consistent with FUV observations. The de- 
tection of a hard X-ray tail from other near-by MSPs such 
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